Phase diagram for unzipping DNA with long-range interactions 
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We present a critique and extension of the mean-field approach to the mechanical pulling tran- 
sition in bound polymer systems. Our model is motivated by the theoretically and experimentally 
important examples of adsorbed polymers and double-stranded DNA, and we focus on the case in 
which quenched disorder in the sequence of monomers is unimportant for the statistical mechanics. 
We show how including excluded volume interactions in the model affects the phase diagram for the 
critical pulling force, and we predict a re-entrancy phase at low temperatures which has not been 
previously discussed. We also consider the case of non-equilibrium pulling, in which the external 
force probes the local, rather than the global structure of the dsDNA or adsorbed polymer. The 
dynamics of the pulling transition in such experiments could illuminate the polymer's loop structure, 
which depends on the nature of excluded volume interactions. 
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Recent experiments using micromanipulation tech- 
niques to pull and stretch single molecules have led the 
way to a better understanding of statistical properties of 
polymer systems that cannot be probed in the bulk jlj. 
Theoretical work on these systems has focused on the 
complex phenomena arising from heterogeneity in the 
sequence of monomers [0. In this Letter we study the 
simpler case of a homopolymer under the influence of an 
external force. We find that some of the same phenom- 
ena caused by disorder in a heteropolymer can arise in a 
homopolymer system near its thermal denaturation tem- 
perature due to inhomogeneities in the polymer's struc- 
ture. We also predict a re-entrant phase in the critical 
force phase diagram at low temperatures, which should 
be observable in both homopolymer and heteropolymer 
systems. 

The Hamiltonian for a polymer in an external field 
takes the general form 0] 
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with d the spatial dimension and b the Kuhn length. 
Throughout this Letter, we set fcs = 1. </>(r) is an at- 
tractive potential arising, for example, from the interac- 
tion with an adsorbing surface, or from the attraction 
between two strands in double-stranded DNA (dsDNA) . 
We assume that <f> is short-ranged, and depends only on 
the local coordinate, r(n), describing the distance from 
the adsorbing surface or the separation between the two 
dsDNA strands. In the case of a homopolymer, 4> does 
not depend on the position along the sequence. Finally, 
the volume interactions between different segments of the 
chain are represented by the term E[r(n)]. 

In the mean-field approximation, interactions between 
different sections of the polymer chain, together with 
whatever external fields are present, contribute to an 
overall effective potential, (f> e e-, for each monomer. This 



field is derived using conditions of self-consistency, and 
for low temperatures it should admit bound states cor- 
responding to adsorbed or zipped conformations. In 
this framework, every monomer far from the ends of 
the molecule experiences the same field, and thus con- 
tributes the same amount to the overall bound-state 
free energy. The mean-field free energy thus reduces to 
J-mf = NX, where A is determined by the self-consistent 
field, o ff 0, ||. Recent theoretical work has explored 
the effect of an external pulling force, represented by a 
vector potential term in the energy, on the mean-field 
thermodynamics [Q. A polymer that is adsorbed to an 
attractive surface or "zipped" together with a comple- 
mentary chain as in dsDNA, can be pulled out of the 
bound state by a force exerted, for example, by a glass 
bead attached at one end of the molecule. The desorbed 
section of a polymer or the unzipped part of dsDNA has 
a force-dependent free energy, T = Ng(F), which de- 
creases monotonically with F. The functional form of 
g(F) can be obtained from a Gaussian or a freely-jointed 
chain model for the stretched polymer. At a critical force 
defined by g(F c ) = A, the adsorbed or double-stranded 
state becomes unstable with respect to the stretched part 
of the molecule, and a sharp unzipping/desporption tran- 
sition occurs . 

The mean-field model relies on the strong assumption 
that every bound monomer experiences the same effec- 
tive field, so that the interaction between the stretched 
and bound parts of the molecule is independent of the 
number of monomers already pulled out. This feature of 
the mean-field model breaks down, for example, in the 
presence of variations in the binding energies for differ- 
ent base pairs in dsDNA, or sequence disorder in an ad- 
sorbed polymer. The effect of such sequence heterogene- 
ity on the pulling transition has been the subject of care- 
ful study, and several interesting phenomena have been 
predicted to emerge in these systems Q. Yet the sim- 
ple mean-field model of the pulling transition can break 
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down even for totally homogeneous polymers when long 
length-scale correlations in the structure appear near the 
temperature of thermal unbinding (denaturation or des- 
orption), T c . Intuitively, thermal fluctuations near T c 
cause large loops and sub-chains to form. The pulling 
force interacts with the bound polymer only in the lo- 
cal region near the end of the stretched-out part of the 
molecule, so the local structure is important in deter- 
mining the behavior of the pulled polymer. Spatial inho- 
mogeneities caused by large loops near T c affect both the 
equilibrium and the non-equilibrium behavior of a bound 
polymer in response to a pulling force. We consider these 
effects in turn. 

The pulling force applied to one end of the polymer 
adds a vector potential term in the overall Hamilto- 
nian 0, giving H = Ho + W pu ii, where 



H pu n = -F ■ r(0) = F / 
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We have assumed that the pulling force, F = Fz, is 
perpendicular to the adsorbing surface. Neglecting the 
effect of the finite size of the molecule, we can drop 
the last term. After tracing over all configurations, 
{r(n)}, in the partition function, we are left with the 
following contribution to the free energy from the pulled- 
away part of the molecule 0): g(F) = — F 2 ^ T . Alterna- 
tively, the freely-jointed chain model, which is more real- 
istic at high pulling forces and low temperatures, gives: 

~ Tsinh(Fb/T) ~ 
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In the mean-field framework, the pulling transition oc- 
curs at g(F c ) = A, or 



F c =y/2d\T/b 2 
F c =(T/b) S mh(F c b/T)e-V T 
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The mean-field approximation should hold good at low 
temperatures, since most monomers will be in the bound 
state and only short loops will form. As T — > 0, the 
free energy per monomer, A, converges to some con- 
stant value, — Ao, corresponding to the binding energy 
of one monomer. The Gaussian model thus predicts 
F c ~ (AoT) 1 / 2 ; however, in this regime the Gaussian 
model for the free energy of a stretched chain is not valid, 
since there is little elasticity in the polymeric bonds. We 
instead rely on the freely-jointed model, which gives at 
low temperatures: g(F) ~ -Fb -Thx{T/2Fb), so that 
F C {T = 0) = A /6. The limiting behavior for T « X Q b is 
given by 



(4) 



We thus expect to observe two desorbed phases, one at 
high temperature and another at sufficiently low temper- 
ature, for F > Xo/b. At high temperatures, the entropic 



advantage of the unbound monomers overrides the ener- 
getic stability of the bound phase. At low temperatures 
the free energy per monomer of the stretched part of the 
polymer or DNA increases as the chain becomes stiffer. 
The effect of lowering temperature near T = is to stiffen 
the chain, making it easier to pull out the polymer or un- 
zip the dsDNA. The temperature at which this re-entrant 
phase should exist is determined by the specific proper- 
ties of the chain, in particular the Kuhn length, b. In fact, 
at low temperatures, we expect F C (T, b) = ^F C (T, b = 1). 
This prediction may allow for experimental determina- 
tion of the relative chain stiffness, related to b, for dif- 
ferent molecules, through careful low-temperature inves- 
tigation of F c . 

Close to the thermal unbinding temperature, T c , as 
A approaches zero, the mean-field approximation breaks 
down because of long-range correlations (corresponding 
to long loops) in the structure of the bound chain. In- 
deed, within the mean-field framework the mean size of 
loops grows near T c as @: (k) ~ A" 1 ~ |T - T c |~ 2 . 
Forming a loop carries an entropic cost, which depends 
non-linearly on the loop size, not simply on the local 
density of monomers. Dealing directly with the interac- 
tion energy, B[r(n)], complicates the analytic problem, 
although attempts have been made to simulate this full 
model numerically ||. We will follow a different route 
first proposed by Poland and Scheraga fft> f|, in which 
the mean-field model is extended by introducing a gen- 
eral form for the entropy of a loop of length k, 



S(k)=A + kln(s)-cln(k), 
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where s is a non-universal constant, and c is determined 
by the properties of the loops. In particular, for non- 
interacting monomers (E[r(n)] = 0) the loops are ran- 
dom walks in d dimensions, with c = d/2, whereas ex- 
cluded volume interactions between monomers tend to 
increase the value of c || . 

The binding energy due to the potential, 0, leads to a 
statistical weight u> = e~^°/ T for each bounded base-pair 
or adsorbed link. In fact, the only relevant parameter in 
this model is the energy difference between the bounded 
and unbounded links, ln(a->) — ln(s), so that we can simply 
put s = 1. 

We work with the grand partition function: Z{z) — 
Sm=i Z(m)z m , where Z(m) is the canonical partition 
function for a chain of m links, and z is a fugacity which 
controls the size of the polymer. In this ensemble, the 
chain can be regarded as an infinite collection of loop- 
train pairs, giving: 



Z(z) = 



Vo(z)U L (z) 
1 - V(z)U(z) 
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where U(z) = £~ x ^z k and V(z) = £~ x 
boundary conditions are Vq = 1 + V(z) and Ul = 1 + 
U(z). We will be interested in the thermodynamic limit, 
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FIG. 1: Phase diagram of the critical pulling force, F c . The 
force is given in arbitrary units, and we have set Ao = 0.4. 
The curves correspond to c = 1.25 (solid), c = 3/2 (dashed) 
and c > 2 (dot-dashed). The area below each curve is the 
bound phase, while the regime above the curve corresponds 
to desorbed or unzipped conformations. 

where the total number of links, (N) — d d ^ z z , diverges. 
From Eq. (|) we see that this limit corresponds to the 
value z = z*, where U(z*) = 1/V(z*) = 1/uiz* - 1. 
The thermal desorption or denaturation transition at T c 
corresponds to z = z c = 1 /s |j| || . 

We can introduce an external pulling force through an 
extra boundary condition: Zp — Vp(z)Z, where Vp(z) — 
1 + Y^kLi^F^ > an( i u f = e~ g ( F '' T is the statistical 
weight per monomer of the stretched part of the polymer. 
The number of monomers pulled away goes like (to) = 
gT^f^j ~ (1 — i-^fz*)^ 1 , so that we naturally identify the 
point wp c = l/z* as the mechanical pulling transition. 
For T close to T c this gives, 

F c ~ (T/b)V-2d\nz* (Gaussian) 
F c ~ (T/b)y/l - z* (Freely- Jointed) 

At the thermal desorption transition, z* — ► z c = 1/s = 1, 
so that F c — > 0. 

The approach to T c is described by the fraction of 
bounded segments (|], 9 = d J" n ^ . We use the fact that, 
near T c and with c > 1, we have \U(z c ) — U(z*)\ ~ 
\z* — z c | , where £ = min(l,c — 1). Thus we obtain 
\z* - z c \ ~ {T-Tcl 1 ^, so that 

F C ~|T-T C |* (8) 

for both the Gaussian and the freely-jointed models. We 
illustrate this result, together with the low-temperature 
behavior from Eq. (|4|), in Fig. |l|. 

Note that at c = 3/2, corresponding to the entropy 
of non-interacting loops, Eq. (||) gives F c ~ \T — T c \, 
in agreement with the mean- field result §, F c ~ VX 
For c > 3/2, F c has a kink at T c , while for c < 3/2, F c 
approaches zero with a continuous derivative. Thus the 
character of the mechanical desorption transition changes 



at c = 3/2, in contrast with the change in the thermal 
desorption transition, with 6 as its order parameter, at 
c = 2. Note that the low-temperature behavior is un- 
affected by the exponent, c, since loop entropy plays no 
role in this regime. 

Although the value of c affects the order of the ther- 
mal denaturation transition J8|, the mechanical desorp- 
tion transition at F c is always first-order in the case of 
a homopolymer. For low pulling forces, F < F Cl there 
is a free-energy advantage to keeping as many links as 
possible bound, while for F > F c the molecule would 
like to unbind as many links as possible. The desorbed 
section of the chain is regarded, in the grand canonical 
ensemble, as a separate system described by the parti- 
tion function, Vf(z*). Thus we recover the mean-field 
results g: (to) ~ |F-F C | _1 , and ((Am) 2 ) - |F-F C |- 2 . 

Our phase diagram, Fig. |l|, agrees qualitatively with 
very recent results for pulled polymers modelled as self- 
avoiding walks on a hyper-cubic lattice in d dimen- 
sions H . Eq. (||) suggests an empirical test for the expo- 
nent, c, based on a precise measurement of the behavior 
of the critical force near T c . In particular, experiments 
on the pulling of DNA and of adsorbed homopolymers 
are expected to display different scaling of F c near T c . 

So far, we have considered the ensemble where the ex- 
ternal parameters, F and T, are fixed, so that the pulled 
polymer remains in equilibrium during the experiment. 
In fact, as long as these parameters are varied sufficiently 
slowly, the system will adjust to changes "adiabatically." 
Thus, for example, the parameter m will be character- 
ized by the force-dependent distribution described above, 
while the globular part of the polymer will display a cor- 
responding, temperature-dependent distribution of loop 
sizes. The behavior of this distribution will depend on c 
close to T c . 

Fast variations in the pulling force, on the other hand, 
may lead to non-equilibrium behavior, in which the dis- 
tribution of loop sizes plays a crucial role. We can es- 
timate the relevant timescale for this "non-adiabatic" 
pulling from the Rouse model for the dynamics of an 
ideal Gaussian chain. The relaxation time for corre- 
lations between the ends of a chain of length N is 
tff = N 2 b 2 t;/(3n 2 T), where £ is the coefficient of friction 
for monomers in the solvent ||. In general, the relax- 
ation time for a section of the chain of length k scales as 
fit ~ k 2 . If the external force dislodges monomers at a 
rate of 1/tF, we see that for tp 3> t/v the polymer will 
remain in equilibrium during the pull-out. 

On the other hand, if tp is much smaller than the 
relaxation time of the smallest loops, the polymer will 
not have time to re-arrange itself as F is ramped up, and 
whole loops will be pulled out at once. In this case, the 
force-extension curve will display jumps corresponding to 
the pullout of individual loops. This behavior is similar 
to that derived in models for polymers with sequence 
disorder B, but the randomness in this case arises from 
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FIG. 2: The critical force, F c , as a function of temperature 
near T c = 1, in the non-adiabatic regime, where the pulling 
timescale is less than the relaxation time of the polymer loops. 
The curves correspond c = 1.25 (solid), c = 3/2 (dashed) and 
c > 2 (dot-dashed). Note that in the non-adiabatic regime, 
F c vanishes below T c . The form of these curves depends on 
the value of Ao, the binding energy of a single link, which is 
here chosen arbitrarily. 



the distribution of loop sizes. The loops serve as a source 
of quenched disorder in this regime, since they do not 
have time to adjust to changes in the mechanical force. 
Alternatively, the distribution of loop sizes could become 
important even for "adiabatic" pulling if the length of 
the stretched part, (m), is less than the typical size of 
loops, (fc). For c < 2, we expect the correlation length, 
and hence (fc), to diverge near T c , so that for fixed F 
the loop size will eventually grow larger than (m), with 
a crossover at (m) sa (fc). 

Thus for either a fast-acting external force, or for T 
sufficiently close to T c , the pulling force interacts with 
the polymer locally, at the level of an individual loop, 
rather than with the entire molecule. Near T c most of the 
monomers are unbound, and we can assume the last loop 
is attached to the surface or the complementary strand 
with only one link. If the loop contains fc unbound links, 
its free energy per link is simply r/k = j n ^ ms 
regime we can write a "semi-microscopic" model for the 
free energy as a function of the number of pulled-out 
mon omcre &: 



Eq. (^|) suggests that the critical pulling force in the 
non-adiabatic regime, F c , is smaller than in the adiabatic 
regime. This is because the configurational entropy of a 
single loop is smaller than that of an ensemble of loops 
with variable lengths. In particular, 
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where the brackets indicate averaging over the distribu- 
tion of loop sizes, (rj k ) oc Y%Li ln "+ s(fc) IhL. Fig. g 
shows the phase diagram in the non-adiabatic regime 
near T c , for three values of c. In all three cases, F c van- 
ishes at a temperature less than T c , in contrast with the 
adiabatic case. This temperature is determined by the 
binding energy, Ao = — Tina;, of the single bond at the 
loop's end. 

Note that, as the mean loop size diverges for c < 2, 
the timescale defining adiabatic pulling diverges. In this 
regime even relatively slow variations in the pulling force 
could be faster than the relaxation time of the loops. 
Simulations and analytical calculations on the equiva- 
lent of Eq. (||) for heteropolymers have shown how the 
size of jumps and plateaux in the force-extension curve 
vary as F approaches F c |Q. As mentioned above, anal- 
ysis of Eq. (||) is complicated by the non-Gaussian dis- 
tribution of the "noise" term, r]k, but simulations could 
elucidate the analogy between the effects of structural 
inhomogeneities in pulled homopolymers and sequence 
disorder in heteropolymer systems. 

The authors would like to thank D. Mukamel and D. 
Lubensky for important comments on an earlier version 
of this paper. 
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where rjk(i) is the free energy per monomer of the loop 
containing link i. The distribution of loop lengths is sim- 
ply P(k) oc e hrik — . Eq. (|9|) resembles the Langevin 
equation recently proposed as a model for the pulling 
transition in a heteropolymer with sequence random- 
ness , but in this case the randomness due to the loop 
distribution has a non-Gaussian distribution, and thus 
cannot be analyzed in the simple framework of Brownian 
motion. 
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